Structural and dielectric properties of Sr 2 Ti04 from first principles 
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We have investigated the structural and dielectric properties of S^TiCU, the first member of the 
Sr n +iTi„03n+i Ruddlesden-Popper series, within density functional theory. Motivated by recent 
work in which thin films of S^TiCU were grown by molecular beam epitaxy (MBE) on SrTiOs sub- 
strates, the in-plane lattice parameter was fixed to the theoretically optimized lattice constant of 
cubic SrTiC>3 (n=oo), while the out-of-plane lattice parameter and the internal structural parame- 
ters were relaxed. The fully relaxed structure was also investigated. Density functional perturbation 
theory was used to calculate the zone-center phonon frequencies, Born effective charges, and the 
electronic dielectric permittivity tensor. A detailed study of the contribution of individual infrared- 
active modes to the static dielectric permittivity tensor was performed. The calculated Raman and 
infrared phonon frequencies were found to be in agreement with experiment where available. Com- 
parisons of the calculated static dielectric permittivity with experiments on both ceramic powders 
and epitaxial thin films are discussed. 

PACS numbers: 77.22.Ch, 63.20.Dj, 61.66. Fn 
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I. INTRODUCTION 

Nonstoichiometric SrTiC>3 resists the formation of 
point defects by forming crystallographic shear (CS) 
phases. 1 For excess SrO, these CS phases form the 
Sr„ + iTi„03„ + i Ruddlesden-Popper (RP) series. The 
structures of the members of the RP series can be viewed 
as a stacking of SrO-terminated SrTiC>3 perovskite [001] 
slabs with relative shifts of (ao/2)[110], the slabs in the 
nth member of the series having a thickness of n cubic 
perovskite lattice constants. This representation of the 
structure for n = 1 is shown in Fig. • 

Considerable fundamental and practical interest in the 
RP series arises from the dielectric properties of its end 
member SrTi03 (n = oo). Strontium titanate is an 
incipient ferrroelectric with a high dielectric constant 
that can be readily tuned by applying a small DC bias. 
SrTiC>3 has low dielectric loss at microwave frequencies 
but a large temperature coefficient of dielectric constant 
(TCF). The prospect of "engineering" the properties of 
the constituent SrTiOs layers by varying n has stimu- 
lated the synthesis and study of the dielectric proper- 
ties of various members of the RP series. Through the 
use of conventional ceramic processing techniques, single 
phase Sr 2 Ti0 4 (n = l), 1 : 2 : 3 ! 4 : 5 Sr 3 Ti 2 7 (n = 2),Wi 5 i 6 



and Si^TisOio (n = 3)j» were formed. Recently, thin 
films of the first five members of the RP series were 
grown by molecular beam epitaxy (MBE) on SrTi03 sub- 
strateSi^S Sr2Ti04, Sr3Ti207, and S^TisOio were found 
to be nearly single phase while SrsTi40i3 and S^TisO^ 
films showed noticeable antiphase boundaries and inter- 
growth defects^ The dielectric properties of the n = 1,2, 
and oo ceramic samples have been characterized by var- 
ious groups 4i£iLi2i The dielectric properties of the n = 3 
and 4 ceramic samples have also been studied but these 
were believed not to be single phased It was found that 
SrTiC>3 had both the highest dielectric constant, e r , and 
the largest TCF of the RP series. Conversely, Sr 2 Ti04 



was found to have the lowest e r of the series but also the 
lowest TCF (an order of magnitude smaller than SrTiOs) 
while having a dielectric loss comparable to SrTiOs. For 
thin films, Haeni et al& performed low frequency and 
microwave experiments to measure the dielectric permit- 
tivity of Sr2TiC>4 and of Sr3Ti207, to be discussed below. 

Here we take the first step in understanding the dielec- 
tric behavior of the Sr„ + iTi„03„+i Ruddlesden-Popper 
(RP) series from first principles by presenting our study 
of the structural and dielectric properties of Sr2Ti04 
[n = 1) using density functional theory (DFT) struc- 
tural optimization as well as density functional pertur- 
bation theory (DFPT). DFT and DFPT have proved to 
be useful tools in investigating the structure, dynami- 
cal and dielectric properties of metals and insulators, in- 
cluding complex oxides (see Ref. for a review). In 
Section El we give the details of our DFT/DFPT cal- 
culations and describe various constraints that we im- 
posed on the Sr2TiC>4 structure. In Section IIIII the re- 
sults of our calculations for the ground-state structural 
parameters, Born effective charges, zone-center phonons 
and electronic and static dielectric permittivity tensors 
at T=0 are presented and compared with experiment. 
In Section IIVI we discuss the sensitivity of the dielec- 
tric response to structural constraints, and through com- 
parison with SrTiOs, explore the key role played by the 
Ti-0 chain geometrjsi^ in the dielectric response and its 
anisotropy in the RP phases. Finally, in Section Ivl we 
summarize our results and main conclusions. 



II. METHOD 

A. First principles calculations 

First-principles density-functional calculations were 
performed within the local density approximation (LDA) 
as implemented in the ABINIT packaged The exchange- 
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correlation energy is evaluated using the Teter rational 
polynomial fit to the Ceperley- Alder electron-gas data, — 
Teter extended norm-conserving pseudopotentials were 
used with Ti(3s 3p 4s 3d), Sr(4s 4p 5s), and 0(2s 2p) 
levels treated as valence states. The electronic wavefunc- 
tions were expanded in plane waves up to a kinetic energy 
cutoff of 35 Ha. Integrals over the Brillouin zone were 
approximated by sums on a 6 x 6 x 2 mesh of special 
/c-points— 



B. Structural Constraints 

We first performed full optimization of the lattice pa- 
rameters and internal coordinates in the reported struc- 
ture of Sr 2 Ti04)£ the body-centered tetragonal K2N1F4 
structure (space group IA/mmm^D^), with the prim- 
itive unit cell containing one formula unit as shown in 
Fig. ^ Ti atoms occupy Wyckoff position (2a) , O x and 
O y atoms (4c), and Sr and O z atoms (4e), the latter 
with one free parameter (displacement along z) each (see 
Table P). This yields the predicted structure under zero 
stress, roughly corresponding to bulk ceramic powders 
and relaxed epitaxial films. The structural relaxation was 
performed using a modified Broyden-Fletcher-Goldfarb- 
Shanno (BFGS) algorithm— to optimize the volume and 
atomic positions, followed by an optimization of the c/a 
ratio. This procedure was repeated to ensure conver- 
gence. The residual Hellmann-Feynman forces were less 
than 2 meV/A. A second structure was considered to in- 
vestigate the effects of epitaxial strain induced by the 
SrTiC>3 substrate on a fully coherent thin film. With 
the same space group, we fixed the in-plane lattice pa- 
rameter to that of SrTiC>3 calculated within the present 
theory (a= 3. 846 A) jil while optimizing the other struc- 
tural parameters: the c lattice constant and displace- 
ments along [001]. This was followed by ionic relax- 
ation until the residual Hellmann-Feynman forces were 
less than 1 meV/A. 

The next section will show that the fully optimized lat- 
tice parameters are underestimated relative to the exper- 
imental values, a common feature in LDA calculations. 18 
To investigate whether our calculations of the dielectric 
response of Sr2Ti04 are sensitive to this underestimation 
of the volume, two additional structures were considered. 
First, with the same space group, we optimized all inter- 
nal coordinates fixing both the a and c lattice constants 
to values obtained by uniformly expanding the lattice 
constant of the in-plane constrained structure using the 
measured thermal expansion coefficient of SrTiOgi^iSli 
(a «7xl0 -6 /K) and a temperature increase of 300 K. 
We also optimized internal coordinates fixing both the a 
and c lattice constants to the values experimentally ob- 
tained for the thin film sample. In both structures the 
ions were relaxed until the residual Hellmann-Feynman 
forces were less than 1 meV/A. 



C. Linear response calculations 

Linear response methods provide an efficient means for 
computing quantities that can be expressed as derivatives 
of the total energy E with respect to a perturbation, 
such as that produced by displacement Ui a of an atom 
or a homogeneous electric field £ . Examples computed 
in this work include the force constant matrix elements 

d 2 E 



dUiaduj/3 

where Ui a is the displacement of atom i in cartesian di- 
rection a from its position in the equililbrium crystal 
structure, the Born effective charge tensor 

8 2 E 





where di a is the uniform displacement of the atomic sub- 
lattice i in cartesian direction a from its position in 
the equilibrium unit cell, with derivatives taken in zero 
macroscopic field, and the electronic susceptibility tensor 



Xaf3 



d 2 E 
d£ a d£/3 



which is related to the electronic dielectric permittivity 
tensor Coo by £oo = 1 + ^X- In this work, we com- 
pute these derivatives using the variational formulation 
of density- functional perturbation theory (DFPT), im- 
plemented in the ABINIT package. 13 i 21 i 22 

The static dielectric permittivity tensor eo can be ob- 
tained directly from the quantities computed by DFPT. 
In general the static dielectric permittivity tensor of a 
nonpolar material can be writtenSiSi as the sum of the 
electronic dielectric permittivity tensor and a sum 
of contributions Ae m from each of the zone-center polar 
modes m: 
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c q/3 — C Q/3 ~ 

where Ae m is given by^Si 







Aire 2 Z*Z. 
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Here, V is the volume of the primitive unit cell and 
Mo is a reference mass taken as 1 amu. MaV ™ a ™0 
can be thought of as an effective plasma frequency, ft 2 ml 
of the m th normal modef^LSS while u> m is the frequency 
of vibration of normal mode m. 



Z^a, which has been 
referred to as a mode effective charge ji^S*^ is given by 



E 



, i' M ° 
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Cm («7) 



(3) 



where £ m is the dynamical matrix eigenvector; the corre- 
sponding real space eigendisplacement of atom i along f3 
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is given by U m (if)) = t\ m /M^ 2 m& Thus we see that 
large lattice contributions to the static dielectric permit- 
tivity tensor are expected if the relevant mode frequencies 
are very low and/or if the effective plasma frequencies are 
large (reflecting large mode effective charges). 

This formalism has been applied in a first-principles 
context in numerous calculations of the zero-temperature 
static dielectric response (for a review, see Rcf. 30) in sim- 
ple and complex oxides. For example, recent calculations 
for zircon (ZrSiOzL^ and zirconia (Zr0 2 )^^ yield good 
agreement with experiment. 

III. RESULTS 

A. Crystal Structure of Sr 2 Ti0 4 

The reported crystal structure of Sr2Ti04 can also be 
viewed as a stacking of TiC>2 and SrO planes along [001] 
as shown in Fig.^b), the stacking sequence being TiC>2- 
SrO-SrO where the second SrO layer is shifted with re- 
spect to the previous SrO layer by ^ao along [110]. In 
previous work (Ref. 35), the internal structural parame- 
ters are presented as v! and v, the distance between Ti02 
and SrO planes and between successive SrO planes, re- 
spectively, and S, the distance along c between Sr and O 
in the same SrO layer, which quantifies the "rumpling" 
of the layer. 

In Table ITTI we present the theoretical lattice parame- 
ters for the structures discussed in Section [H] and com- 
pare them with the experimental values for ceramic sam- 
ple a 2 ' 4 '^ 52 and for thin films epitaxially grown on SrTi03 
(a=3.905A) substratesi&& The lattice constants of the 
fully relaxed structure, Table ITlTa'). are smaller than the 
measured lattice constants of the ceramic powder. Specif- 
ically, the a lattice parameter was calculated to be less 
than experiment by 1.5%, which is within the error typ- 
ically associated with the LDA, while the c lattice pa- 
rameter is underestimated by 2.2%. The smaller value of 
c measured for the thin film samples reduces this dis- 
crepancy to 1.1%. There is considerable rumpling of 
the SrO layers, Ss r o = 0.05ao- This rumpling is such 
that the O z and Sr atoms move in a direction away from 
and towards the Ti02 layers, resulting in a Ti-O z bond 
slightly larger (2.9%) than the Ti-O y (or equivalently Ti- 
O x ) bond. As pointed out by Noguera, Ruddlesden and 
Popper had assumed that this rumpling was equal to 
zero. The present study supports the previous calcula- 
tionS by finding a nonzero rumpling. The experimental 
work of Venkateswaran et alS^ also suggested a nonzero 
rumpling of the SrO layer (t>s,.o=0.10)r& but found that 
the Ti-O z bond is shorter than the Ti-O y bond. 

Next, we consider the structure in which we con- 
strained the in-plane lattice constant to that of theoret- 
ical SrTi0 3 (a = 3.846 A), Table E3b). This places the 
system under tensile in-plane stress equal to that of an 
experimental sample of Sr2Ti04 coherently matched to 
SrTiOa (lattice mismatch 0.6%). As can be seen in Ta- 



ble [D] this has almost no effect on the optimized values 
of the other structural parameters, including c. 

Finally, we consider the optimized structural param- 
eters for the two expanded structures. The slight ther- 
mal expansion of the in-plane constrained structure has 
a correspondingly slight effect on the internal parame- 
ters, while the expanded c of the experimental thin film 
structure has a larger effect. We see that a roughly ho- 
mogeneous expansion of the lattice, i.e. column (a)— »(d) 
and (b)— >(c), decreases the rumpling of the SrO layer. 



B. Born Effective Charge Tensors 

In the Sr2Ti04 structure, the site symmetries of the 
Ti atom, occupying Wyckoff position 2a, and the O z 
and Sr atoms, both occupying Wyckoff position 4e, are 
tetragonal, while that of the O x and O y atoms. Wyckoff 
position 4c, is orthorhombic. As a result, all Z*'s are 
diagonal with two (Z*(Ti), Z*(O z ), and Z*(Sr)) and 
three (Z*(O y ) and Z*(O x )) independent components. 
Table IIIII displays the diagonal components of the cal- 
culated Born effective charge tensors in Si - 2Ti04 for the 
four structures considered. We see that the Z*'s are 
relatively insensitive to the various volume and strain 
constraints imposed, consistent with what was shown 
for isotropic volume changes in BaTiOg^ and for pure 
tetragonal strain in KNbOa4^ Also, we note an anoma- 
lously large value for Z* x (Ti) = Z* y (Ti) (nominal charge 
+4) and for Z* x (O x ) = Z* y (O y ) (nominal charge -2), 
these being the components of the Z*'s corresponding 
to motion parallel to the Ti-0 bond along a direction in 
which the infinite Ti-0 chains have been preserved. In 
contrast, the anomalous parts of Z* z (Ti) and Z* z (O z ), 
which again correspond to motion parallel to the Ti-0 
bond, but along [001] where the Ti-0 bonds do not form 
continuous chains, are found to be less than half of those 
along the continuous Ti-0 chain. 



C. Phonon Frequencies at V 

For the Si~2Ti04 structure, group-theoretical analysis 
predicts that the 18 zone-center optic modes transform 
according to the following irreducible representations 

1 optic 

= 2A lg © 2E g © 3A 2u © AE U © B 2u (4) 

of which the Ai ff and E g modes are Raman active, the 
A2u and E u modes are infrared active, and the B2« mode 
is neither Raman nor infrared active. By using projection 
operator methods, it can be shown that the Ai g , A 2u , and 
B2u modes involve motion along [001] while in the E s and 
E n modes, atoms move along [100] and [010]. A complete 
listing of the symmetry adapted lattice functions is given 
inRefs.liJandlii 

Table llVl displays our calculated frequencies and mode 
assignments for the four structures considered. Compar- 
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ing the calculated phonon frequencies for the fully re- 
laxed structure (a) with those measured by Burns et al. 
on ceramic samples, also given in Table HH we find ex- 
cellent agreement. The calculations allow us to confirm 
the mode assignments of the A 2u (TO\) and E U (T03): 
Burns' suggestion of assigning the lower mode at 242 
cm" 1 symmetry and the higher mode at 259 cm -1 
E„ symmetry appears to be correct. In addition, Burns' 
assessment that the strong LO feature at «440 cm -1 in 
the measured reflectivity spectra was "probably" due to 
both a A 2u (LO) and a E u (LO) mode is consistent with 
our calculations, where we found A 2u {LO) —479 cm -1 
and E„(LO) =451 cm -1 m& In fact an average of these 
two frequencies is within 2 cm" 1 of the quoted value 467 
cm -1 . Finally, we obtain the frequencies of four modes 
that could not be separately identified in the experiment, 
namely the A 2u {T02) mode at 378 cm -1 , A 2u (LO) mode 
at 252 cm" 1 , the E U (T04) mode at 611 cm" 1 , and the 
B 2u mode at 303 cm" 1 . 

In Table El we give the normalized dynamical matrix 
eigenvectors for the fully relaxed structure (a). We see 
that the lowest frequency A 2u (T01) and E U (T01) modes 
involve Sr atoms moving against a fairly rigid TiC>6 octa- 
hedron, with larger deformation of the octahedron for 
the E M mode. The highest frequency A2 U (T03) and 
E u (T04) modes each involve large distortions of the oxy- 
gen octahedra with the apical oxygens moving against 
the planar oxygens in the A 2u mode and a corresponding 
distortion for the E M mode. The calculated A 2u (T02) 
mode involves Ti atoms moving against a nearly rigid 
oxygen octahedron, giving rise to a relatively large mode 
effective charge. Finally, the E„(T02) mode, character- 
ized by Burns et al. as a weak infrared active mode, 
has displacements of O y and O z very close to those of 
the triply degenerate silent mode in the cubic perovskite. 
These patterns are substantially different from those ob- 
tained in model calculations of the lattice dynamics of 
K2ZnF4^i used by Burns to characterize the eigenvec- 
tors of S^TiCU, especially for the highest frequency 
A 2u (T03) and E U (T04) modes and the A 2u (T02) mode. 

The effective plasma frequencies (Sec. Ill CJI are re- 
ported in Table IVII There is relatively little variation 
among the different modes of the same symmetry in the 
same structure, with il p . m «400 cm" 1 and 800 cm" 1 for 
the A2« and E„ respectively, and even less variation with 
changes in the structure. The largest calculated effective 
plasma frequency is obtained for the A 2 „(T02) mode. 
While the corresponding oscillator strength, which goes 
like £lp m , should make this the most prominent mode 
in a single-crystal infrared study, it was not separately 
identified in the ceramic sample of Burns. In ceramics, 
while the reststrahlen bands do not overlap for modes 
transforming as the same irreducible representation, the 
A2 U modes can overlap the E u and vice versa. In fact we 
see that not only the A2 U (T02) mode, but also a second 
high-oscillator strength mode, E tt (T04), not identified in 
the measurements of Burns, lie in the middle of a wide 
reststrahlen band of the other symmetry type. As noted 



by Burns, this would undoubtedly complicate the inter- 
pretation of the reflectivity spectrum and may explain 
why neither of these modes were identified. 



D. Dielectric Permittivity Tensors 

Here we present our calculation of the static dielectric 
permittivity tensor eo for the structures considered in Ta- 
ble [H] By the symmetry of the S^TiCU structure we see 
that the tensors are diagonal and have two independent 
components, along directions parallel to and perpendicu- 
lar to the TiC>2 layers. In Table lVni we display the results 
of our calculations for these three tensors (eoo, (-ionic, and 
eo) for Sr2Ti04 under the various structural constraints 
previously discussed (see the caption of Table [H]). ^ can 
be seen that the static dielectric permittivity tensor is 
quite anisotropic, with the in-plane components nearly 3 
times as large as the component along [001]. 

To improve our understanding of this anisotropy we ex- 
amine the contribution to ti on ic from individual phonon 
modes. In Table IVII we show the effective plasma fre- 
quency, r2 p , m , and the contribution to the static dielectric 

permittivity for each IR-active phonon mode, equal to 

n 2 

J2" 1 ■ The A2„ and the E u modes contribute to the com- 
ponents of eo along directions perpendicular and parallel 
to the TiC>2 layers, respectively. We see that the domi- 
nant mode contributing to the rather large anisotropy is 
the E„(T01) mode, while the largest Q_p, m is associated 
with the A2„(T02) mode. In fact the effective plasma 
frequency of the A 2u (T02) mode is ~40% larger than 
that of the E U (T01) mode. The fact that the frequency 
of this A2j t (T02) mode is more than twice that of the soft 
E n (T01) mode explains its relatively small contribution 
to the dielectric permittivity tensor. 

Finally, to compare the calculated dielectric permittiv- 
ity tensors with the measured dielectric constant of the 
ceramic samples, e r , we compute an orientational aver- 
age, (average: of our calculated results. From Table IVTTl 
we see that (average = 38 for the fully-relaxed structure 
agrees quite well with e r = 34 — 38 measured by the var- 
ious groups. The comparison with the thin film results 
is more problematic and will be discussed in the next 
section. 



IV. DISCUSSION 

In this section, we discuss the main features of the T=0 
dielectric response of Sr2Ti04 and its relationship to that 
of SrTiC>3. We show that the anisotropy and sensitiv- 
ity to changes in the lattice constants through epitaxial 
strain and thermal expansion can be understood by fo- 
cusing on the Ti-0 chains and chain fragments. Finally, 
we discuss the comparison of the results of our calcula- 
tions to available experimental data. 

As noted in the previous section, the static dielectric 
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tensor of Sr2Ti04 is rather anisotropic, with the response 
in the planes of the layers about three times greater than 
the response perpendicular to the layers. This can be 
understood by analyzing the Ti-0 chains. In SrTiOs, 
the large dielectric response can be directly linked to the 
presence of infinite Ti-0 chains running in all three carte- 
sian directions, with anomalously large Born effective 
charges and low frequencies for modes with Ti displac- 
ing relative to O along the chain 55 ; 56 (see Table IVIIIl! . 
In Sr 2 Ti04 the infinite Ti-0 chains lying in the TiC>2 
planes are preserved, while the Ti-0 chains along [001] 
are broken into short O-Ti-0 segments by the relative 
shift of the SrO-terminated perovskite slabs (see Fig- 
ure E^c)). Correspondingly, as can be seen in Table ITTTl 
the Z*'s in Sr2Ti04 for Ti and O displacing along the 
infinite chains in the Ti02 layers, Z* y (Ti) (Z* x (Ti)) and 
Z* y (O y ) (Z* x (O x j), are only slightly smaller than those 
of SrTiOs. 45 In contrast, the breaking of the Ti-0 chains 
into O-Ti-0 segments along [001] has a dramatic effect 
on the anomalous component of the Ti and O Born effec- 
tive charges for motion along the segment, Z zz (Ti) and 
Z* zz (O z ), reducing them by over a factor of 2. While it 
is true that the local anisotropy caused by the somewhat 
larger Ti-O z bond length compared with the Ti-O y bond 
length will also reduce the corresponding Born effective 
charges, this is a much weaker effect as suggested by the 
following calculation. We fixed the Ti-O z distance to that 
of the Ti-O y bond length of structure (b) (where the in- 
plane plane lattice parameter was fixed to the lattice con- 
stant of SrTiOs), creating a structure whereby the Ti-0 
distances are equal in all three cartesian directions. The 
Born effective charges along [001] were found to increase, 
but only slightly (Z* z (Ti)=5.29 and Z z * 2 (O 2 )=-3.80). 

This reduction of the anomalous component of the 
Born effective charges is quite similar to the reduction 
of the Z*'s observed in BaTiOa due to the ferroelec- 
tric transition. Indeed, Ghosez et alM^ found that 
Z*(Ti) = 7.29 and Z*{0\\) = -5.75 in cubic BaTi0 3 
while in the tetragonal phase they found Z* y (Ti) = 6.94, 
Z* y (O y ) = -5.53, Z* zz (Ti) = 5.81 and Z* zz (O z ) = -4.73. 
It was explained that this reduction of the Z*'s going 
from the cubic to the tetragonal phase resulted from the 
displacement of the Ti atom along the ferrroelectric axis, 
resulting in a series of long-short Ti-0 bonds thereby 
"breaking" the Ti-0 chains^ In Sr 2 Ti04 this breaking 
of the Ti-0 chains can be thought of as being caused not 
by alternating Ti-0 bond lengths along a given direction, 
but by the shift of the SrO-terminated SrTiOs slabs, re- 
sulting in a SrO-SrO anti-phase boundary perpendicular 
to [001] as discussed in Sec. 1111 Al Finally, this breaking 
of the chains in Sr2Ti04 is arguably a much "stronger" 
effect than the cubic to tetragonal transition in BaTiOa, 
as evidenced by the greater reduction of the Born effec- 
tive charges along [001]. 

The breaking of the [001] Ti-0 chains into segments 
lowers the dielectric response relative to SrTi0 3 also by 
stiffening the low-frequency mode associated with rela- 
tive displacement of Ti and O along the chains. The 



in-plane E u (TOl) mode in Sr2Ti04 for the fully relaxed 
structure (a) is at 148 cm -1 (134 cm -1 for the structure 
with in-plane lattice constant fixed to SrTiOs (b)) com- 
parable to our calculated value of 103 cm" 1 for SrTiOs, 
while the corresponding out-of-plane A^ u (T02) mode 
is at 378 cm -1 for structure (a) (391 cm" 1 for struc- 
ture (b)). As we discussed in the previous section, the 
large anisotropy is almost completely accounted for by 
the large Ae m of the in-plane E u (TOl) mode as shown 
in Table \VJ\ 

The topology of the Ti-0 chains thus appears to be 
much more important than the bond lengths and the de- 
tails of the geometry. Indeed, our calculations show that 
the dielectric permittivity tensor of Sr2Ti04 is not very 
sensitive to an in-plane tensile strain. From Table IVIII 
we see that e average for the fully-relaxed structure (a) 
and for the structure with the in-plane lattice constant 
constrained to the theoretical SrTiOs value (b) only dif- 
fer by ~15%. Also, we see that this difference originates 
from the component of £i on ic along in-plane directions 
while the lattice contribution to the dielectric permit- 
tivity tensor along [001] remains relatively constant. In 
particular, we see that e^nic along directions in which the 
Ti-0 chains were preserved is relatively more sensitive to 
a fixed strain than along directions in which these chains 
were broken. In fact this sensitivity can be attributed 
to a single E U (T01) mode where it can be seen that 
the sum of the contributions to e^mc from the remain- 
ing E„(TO) modes remains nearly constant for all the 
structures considered. As was pointed out in Sec. IIII Dl 
this E„(T01) mode also contributes predominately to 
the rather large anisotropy of the static dielectric permit- 
tivity tensor. Thus, the difference in response between 
fully coherent and fully relaxed films is expected to be 
relatively insignificant, unlike the case of epitaxial thin 
films of SrTiOs, in which the induced strain would in- 
fluence the lattice dynamics and presumably could have 
significant effect on the dielectric properties' 6 ^ In addi- 
tion, this insensitivity should reduce the TCF, especially 
for e 33 . 

Anisotropy of the dielectric tensor has not been pre- 
viously experimentally determined in Sr2Ti04, though 
measurements have been carried out for single crystals 
of some isostructural oxides CaYA104, CaNdA104, and 
LaSrA104^ In the latter case, the anisotropy is very 
small (and in the case of LaSrA104, the out-of-plane 
component of the dielectric tensor, £33= 20.02, is slightly 
greater than the in-plane component en= e22=16.81). 
The absence of a large in-plane response in these mate- 
rials highlights the unique physics of the Ti-0 chain. 

The anisotropy of the dielectric tensor is an impor- 
tant factor to be taken into consideration when compar- 
ing the computed dielectric response with experiment, 
and indeed in correctly interpreting experimental mea- 
surements. For ceramics, anisotropy of the single crys- 
tal response can give a range of values for the mea- 
sured dielectric response, depending on the size, shape 
and interaction of the grains. In the single crystal films, 
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the anisotropy complicates the interpretation of the mi- 
crowave microscope method of determining the dielectric 
response. The effective dielectric constant obtained by 
this technique is in fact an average of the dielectric tensor 
components determined by the field distribution near the 
probe, which is itself determined by the anisotropic di- 
electric response. While an exact value requires detailed 
modeling of the field configuration, the expectation is 
that the measured value should be an average weighted in 
favor of the high dielectric response components. Thus, 
the value e e ff — 44 obtained by Haeni^ is not £33, but 
is a weighted average of the three components, in accor- 
dance with our calculations. The agreement of this e e // 
with the low-frequency electronic measurement is appar- 
ently coincidental. The latter measurement determines 
633, but at low frequencies it is not unusual for additional 
extrinsic contributions to raise the value of the dielectric 
response. 



available, about the anisotropy of the response. This 
anisotropy is key to interpreting the available experi- 
mental data on the dielectric response of ceramic and 
thin film samples. We identify a dominant role played by 
Ti-0 chains in the lattice contribution to the response, 
which allows us to understand the anisotropy. This pic- 
ture should generalize to higher RP phases of Sr-Ti-O, 
and further could be useful in tailoring the dielectric re- 
sponse of other perovskite-based systems. 
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FIG. 1: (Color online) The structure of Sr2Ti04 (space group IA/mmm) can be viewed as (a) a stacking of SrO-terminated 
SrTiOs perovskite [001] slabs, (b) a stacking of TiCh and of SrO planes along [001], and (c) a series of Ti-O chains, infinitely 
long in the plane along [100] and [010], and of finite extent along [001] (the Sr-atoms have been removed for clarity). 



TABLE I: Wyckoff positions for Sr 2 Ti0 4 . 



Atom Wyckoff position and point symmetry Coordinates 

Ti (2a) 4/mmm 0,0,0 

O x , O y (4c) mmm ±,0,0; 0,±,0 

O z (4e) 4mm 0,0,±z Oz 

Sr (4e) 4mm 0,0,±zsv 
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TABLE II: Structural parameters of S^TiCU. (a) Fully relaxed cell, (b) Relaxation along c-axis, with a lattice constant fixed at 
theoretical SrTi03. (c) Relaxation with lattice constants fixed at values corrected for thermal expansion, (d) Lattice constants 
fixed at experimental (Ref. 8). Lattice constants, bond lengths, and internal parameters are in A. 
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TABLE III: Nonzero components of the calculated Born effective charge tensors Sr2Ti04- 
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TABLE IV: Phonon frequencies (cm -1 ) and mode assign- 
ments for Sr2Ti04. Symmetry labels follow the convention 
of Ref. IB4I Experimental values are from Ref. l4ll 
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TABLE V: Dynamical matrix eigenvectors, £ m , of Sr2Ti04 for the fully relaxed structure with a=3.822A and c=12.32A. The 
corresponding eigendisplacement in real space can be obtained by dividing each value by the appropriate mass factor \fMl. 
Modes Ai g , A 2u , and B 2u involve motion along [001], while the two-fold degenerate modes E g and E u involve motion along 
[100] and equivalently along [010]. 
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TABLE VI: Effective plasma frequency, Q p>m (cm 1 ), and 
mode contribution to the static dielectric tensor, Ae 
Sr 2 Ti0 4 . 

(a) (b) (c) (d) 



Mode 
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Ae m 




Ae m 
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TABLE VII: Nonzero components of the electronic and the ionic dielectric tensors for Sr 2 TiC>4, (e xx , e yy , e zz ). Experimental 
measurements on thin films (Ref. 0) and on ceramic samples (Ref. were conducted at room temperature except for f 

(T=15K.) 





Experimental 










Ref. 9 


Ref. 7 Ref. 4 


Ref. j5 (a) 


(b) 


(c) 


(d) 






( 5.09 5.09 4.81) 


( 5.08 5.08 4.82) 


( 5.07 5.07 4.82) 


( 5.08 5.08 4.79) 






(42.8 42.8 12.5) 


(52.6 52.6 11.9) 


(57.3 57.3 12.2) 


(74.7 74.7 14.2) 


eo 




(48 48 17) 


(58 58 17) 


(62 62 17) 


(80 80 19) 


^-average 44 zb 4 


38 34, 37 1 " 


37 38 


44 


47 


60 



TABLE VIII: The structural and dielectric properties of 
SrTiC>3 calculated within present theory. (6x6x6 grid, 45 Ha.) 

Lattice Constant Phonons (cm -1 ) uj m Q, p , m 

a= 3.846 A A 2u (T01) 103 1275 

Born Effective Charge A 2u (T02) 189 827 

Z*(Tz)= 7.26 A 2u (T03) 587 946 

Z*{Sr)= 2.55 A 2u (L01) 166 

Z*(O x )=-2.04 A 2 „(L02) 451 

Z*(Oy)=-5.72 A 2u (L03) 822 
Dielectric Constant 
e°° = 6.2, e = 181 

Phonon Eigenvector { un, u$ r , Uo\\ > u o±, u o±} 

A 2u {T01) {-0.51, -0.14, 0.37, 0.54, 0.54} 



